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Radio-sonde is used for measuring atmospheric parameters such as temperature, 
relative humdity, and winds. It is essentially an instrument package (sonde) that has 
sensors to measure temperature, humidity, and winds and are attached to a balloon filled 
with lighter-than-air gas such as helium Hence as the balloon goes up in the atmosphere 
it collects data and using a radio it relays the data back to the earth based receiver. We 
have obtained about 2300 such radio-soundings from the TRMM field campaigns held in 
the year 1999. The humidity data obtained from the Vaisala radio-sondes were found to 
be corrupt i.e. during the sounding the atmosphere is actually moist but the radio-sonde is 
measuring it as "dry". In order to remove this "dry-ness" in the humidity data we have 
applied some humidity correction procedure. These procedures were developed at the 
National Center for Atmospheric Research (NCAR) in Boulder, CO based on some 
laboratory tests. Some of the correction procedure involves correcting the data based on 
age of the sonde i.e. error calculated based on the number of days elapsed from the date 
the package was manufactured from the day of its actual use in the field. We have 
interestingly found that sensor arm heating error - one of the problems that was 
independent of age in the correction method of NCAR is actually dependent on the age of 
the sonde. We have found that after correction of the relative humidity data the available 
energy due to convection (the energy that drives the atmosphere i.e. helps cloud 
formation, etc. and as calculated using the radio-sonde data) has increased from before. 



ABSTRACT 


Hundreds of Vaisala sondes with a RS80-H Humicap thin-film capacitor 
humidity sensor were launched during the Tropical Rainfall Measuring Mission 
(TRMM) field campaigns in Large Scale Biosphere- Atmosphere held in Brazil 
(LBA) and in Kwajalein experiment (KWAJEX) held in the Republic of Marshall 
Islands. Using Six humidity error correction algorithms by Wang et al. (2002), 
these sondes were corrected for significant dry bias in the RS80-H data. It is 
further shown that sonde surface temperature error must be corrected for a 
better representation of the relative humidity. This error becomes prominent due 
to sensor arm-heating in the first 50-s data. 
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1. INTRODUCTION 


During the TRMM field campaign, a combined total of 709 soundings 
were taken over the KWAJRX (R/V Ron Brown) and LBA sites of Rebio Jaru 
and Rolim de Moura (see Figures la, and lb). These soundings have undergone 
quality control (Roy and Halverson, 2002) for general use by the scientific 
community. We expect that the corrected sounding s will provide better 
estimate of the latent heating budgets in the tropics especially when compared 
with model-derived and satellite estimates, as well as the characteristics of 
convective systems (see Halverson et al. 2001). During our quality control we 
observed a large dry-bias in the Vaisala soundings (see Figure 2) when 
compared with coincident VIZ soundings launched from nearby stations. To 
correct these biases we used the empirical polynomial fits and regression 
equations based on National Center for Atmospheric Research - Atmospheric 
Technology Division (NCAR/ ATD) series of laboratory tests as reported in 
Wang et al. (2002a), hereafter referred to as W02. We found a consistently good 
response on application of the algorithm to the present set of soundings. 
Corrections are done by applying the equations for the temperature-dependence 
model, modeled ground-check (GC,) contamination, basic calibration and 
daytime sensor arm-heating error (SAHE) . In Section 2 this note we describe the 
sensor RS80-H and how the data were obtained during the TRMM field 
campaigns. Section 3 describes the procedure adopted for correcting the Vaisala 
RH dry bias and the assumptions that we made. Generally it was thought that 
the sensor arm heating biases are independent of the sonde age but after 
application of the correction procedure we have found an inherent age- 
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dependence. A brief description of the results obtained after correction appear in 
Section 4. The purpose of this note is to evaluate the RS80-H correction 
algorithm after application to the TRMM soundings and also to report an age- 
dependent variation of the SAHE which is used in W02 based on Cole and 
Miller (1995). These evaluation results may be useful for development of a 
reference radiosonde system as mentioned in Wang (2002b). Because an overall 
improvement in the RH data is obtained using W02 algorithm hence a brief 
discussion about the resulting data quality based on the site-mean convective- 
available potential energy (CAPE) is also given in section 4. In section 5 we 
briefly discuss about the RH corrected version 'e' Vaisala sounding products 
that have been delivered to the Goddard Distributed Active Archive Center (G- 
DAAC). 

2. SENSOR AND DATA DESCRIPTION 

During the TRMM field campaigns relative humidity was measured by use of 
thin-film capacitive hygristor.. The hygristor measures in the range 0 to 100% 

RH with a 1% resolution. This sensor has a lag time of 1-s with 6 m/ s flow at 
1000 mb and at +20 degree C temperature. Data were obtained using the 
DigiCORA II model MW-15 ground system and processor. GC parameters are 
typically measured by an independent set of sensors at the surface. The GC 
parameters are not used to adjust the sounding's factory determined calibration 
coefficients. For the case of R/ V Ron Brown launches an independent integrated 
suite for Improved Meteorological Instruments (IMET) as described in Hosom et 
al. (1995) was used at the surface for obtaining 1-min average data at sonde 
launch time. Calibration was done by feeding in punched-hole strip during the 
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GC procedure. The GC parameters (pressure, temperature, RH, wind-speed and 
wind direction) are then entered and offsets are computed and recorded in the 
sounding data file. For the LBA launches at Rebio Jam and Rohm de Moura hand 
held psychrometer, thermometer, and wind-measuring instrument was used for 
GC purpose. The Vaisala RH sensor samples approximately every 1.5 sec. A 
least-squares technique is then used to generate smoothed data at 10-s temporal 
resolution. This is done by the Vaisala soimding processing algorithm itself. 

There are in all 297 soundings from R/V RHB in the period between July 25, 
1999 and September 11, 1999. The sondes were launched every three hours 
continuously from R/V RHB but were disrupted in the period Aug 19-27, 1999. 
All R/V RHB data were re-formatted to a consistent LBA format (Roy & 
Halverson, 2002). There were 412 soundings performed from the LBA sites of 
Rebio Jam and Rolim de Moura. Full time operations at Rebio Jam and Rolim de 
Moura did not begin until Jan 24, 1999. Operations were terminated at Rolim de 
Moura on February 21, 1999 since the equipment was moved to another site. 
Operations were terminated at Rebio Jam on Feb 25, 1999. The raw data were 
subjected to visual check procedure and the Skew-T InP charts were generated 
for each sounding and were visually inspected to find errors in the data such as 
super-adiabatic lapse rates, low-level temperature and humidity problems and 
few other problems such as rapid drying, and rapid heating or cooling. 
Preliminary quality flagging of the data sets were performed and version 'd' files 
were created using a method similar to the TOGA-COARE quality flagging 
algorithm as mentioned in Loehrer et al. (1996). The NCAR/ ATD algorithm 
mentioned in W02 for RS80-H correction was applied to these version 'd' files. 
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Missing data have been flagged and no attempt has been made to substitute 
missing data with any other estimate. 

3. METHOD 

An overview of the correction algorithm is given in Fig. 10 of W02. A 
simplified flow-diagram is shown here in Figure 3. A GC correction was 
performed using Equation 4-10 given in W02: 

U' = U„-^U'gc, ( 1 ) 

where u„ is the sonde-measured RH and aUgc is the GC error measured from 
an independent instrument. Wherever the GC error values were not recorded 
we used the modeled GC values given by the sensor aging error (SAE); the SAE 
were determined using Eq. 4-11 of W02 (Eq. 3 below). The resulting GC 
correction is added to the entire sounding RH data and then the ambient 
temperature-dependent error model -1 (TDl), basic-calibration model error (BC), 
chemical contamination correction (CC), temperature-dependent correction 

model -2 (TD2) and 10-s Daytime SAHE, and DSAHlO-s corrections are done 
using the Equations 5-lH, 4-6H, 4-lH, 5-2H, 4-7 and 4-8 respectively of W02 


TD\ = -0.61 + 0.956117' + (0.03 1 + 0.00359U')/ 
+(-0.00033 - 0.0000727 U')t^ 

+[(0.0000014 - (9.6£ - 8)U']P 
+[(-3. IE - 09) + (5.43 1£ - 09)U']t* 

(2) 

5A£ = -0. 1638 + 1 .4766d - 0.2257 

(3) 

BC = Hq + HjU + H2U'^ + 

where 

Ho =-0.3019,//] =-0.0081, //2 =0.0011, an^///3 =-1.23£-05 

(4) 
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CC = [khQ + (khx )d + {kh2)(d^)+ikh^){d^)+(.khi){d*)]x 

[pho + (pfti )(t/) + iph2 )(.U^) + {ph2 )iU^)] 
where, 

kho =0.018866, )t/zi =1.978206, ;t/i2 =-1.34278,^/13 =0.36915724 

A:A4 =-0.03241, /j/iQ =1.6994, pA, = 0. 1 368, ;iA2 =-00018.^3 =1 41 05£- 05 

_ 1/^+0 61-0 031^ + 0.000331^ +0.00000141^ + (3. IE- 09)i'* 

0.960 1 + 0.0003591 - 0.000085/ ^ + (9.3E - 08)/^ + (6.93 1 £ - 09)r " 


(5) 

( 6 ) 


es(.Tss) 


(7a) 


DSAHj = {SAHE)e~’> ' (7b) 

where the GC value of RH ( u’) is used from (1) to compute TDl in (2). The 
ambient temperature / is used from the sounding. The serial numbers for each 
sounding was obtained from the manufacturer (Vaisala, Inc.) . Age (d) of the 
sonde (in year) was calculated using the sonde serial numbers. In this study the 
age of sonde ranged from 0.06 year to 2.07 years. Because the basic calibration of 
each sensor is done during the production process and using an averaged 
calibration model, Wang et al. (2002a) used an accurate test chamber and found 
there were errors associated with the average model at various RH points. 

Hence BC error computed using (4) is used for RH correction which in turn 
become input to the contamination correction model (see Fig. 3). The 
contamination is caused mainly due to occupation of non-water molecules at the 
binding sites hence reducing the ability of water molecule absorption by the 
polymer. This contamination process become prominent with age of the sonde 
due to a longer exposure to the packaging material that out-gasses on to the 
substrate. A polynomial fit for an average approximation of the RH error due to 
this contamination (CC) process with sonde age (d) is mentioned in W02 and is 
given here in (5). After this the CC quantity is brought back to an RH value at 
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ambient temperature using (6) and then the TD and SAE are added to get an RH 
value which is not yet corrected for the radiation induced SAKE. This arm- 
heating error is calculated only for the first 5 levels of sounding. This is limited to 
the first 50-s since the thermal time constant (t) of the RH sensor is 13 seconds 
(as in 7b) and it is expected that the RS80-H will equilibrate within a period 
which is 4 times the thermal constant. SAHE is computed using the surface - 
reference RH (U^r), and using the saturated vapor pressure once calculated using 
the sonde surface level temperature (Tgs) and another time using the surface 
reference temperature measured using an independent instrument (Tgr). The 
saturation vapor-pressure data appearing in (5) are computed using method as 
mentioned in Bolton (1980). The surface RH value as measured by the sonde is 
corrected using the GC procedure and is used to compute the surface vapor 
pressure appearing in the denominator of term 2 of (5). No model-based 
correction is imparted to the sonde surface RH value. It is also assumed that 
there were no occurrence of low-level inversion during the day-time sounding 
periods and since the sonde RH at 10-s and above are corrected using the W02 
algorithm, it is then expected that the error in a well-mixed boundary layer 
(during day time) will only be caused due to radiation-induced sensor arm- 
heating. The choice for use of the sonde registered surface RH value is made 
since the sonde measured surface RH value should ideally remain constant with 
height within the well-mixed boundary layer. For SAHE calculations that are 

separated by approximately 10-s (tj) we have the value of damping factor exp(- 
t;/13) in (7b) ranging from 0.46 to 0.02. 
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4. RESULTS 


a. Overall RH correction 

We have incorporated humidity correction for almost all levels and for all 
soundings obtained from the three TRMM field sites following the algorithm 
W02 as shown in Figure 3. The average daytime, and night-time vertical profiles 
of the uncorrected minus corrected RH values are compared. These profiles are 
shown for each of the sites in Figures 4a, 4b, and 4c. The day-time corrections 
include the SAKE modification to the first 50-s of the sounding, or up to about 
250 m above ground level. From Figures 4 it is evident that moisture loading 
due to the current algorithm is prominent in the mid-levels and reach maximum 
in the 200-100 mb region for all the three sites. Figures 5a and 5b show the total 
counts of observations for various age of the respective sonde. For Rebio Jam site 
maximum observations were made with sondes in the age group 0.08 to 0.12 
years, and for Rolim de Moura site maximum observations were also made in the 
sonde age ranging between 0.08 to 0.12 years but there were few soundings that 
were manufactured almost 2 years prior to their use. Maximum observations 
made from R/ V RHB were obtained from sondes that were manufactured in the 
range 0.19 to 0.33 years prior to their use. 

b. Variation of mean correction parameters with sonde age 

Figure 6 shows the variation of GC, CC, DSAH-lOs, with sonde age, and the 
mean of the uncorrected minus corrected RH as a function of sonde age. Sixth 
order polynomial fit for each type of error is shown in Figure 6. We have the 
following results: (a) The mean modeled GC increases with age, (b) The mean of 
CC correction varies with age, (c) The mean contamination, and mean DSAH- 
10s are found to vary in a cyclic manner with age. The CC correction is 
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dependent on age of the sonde as shown in W02 and Eq. (5) here. However we 
have investigated into the reason why the SAHE and hence the DSAH-lOs vary 
in a cyclic fashion with the age of the sonde. 
c. Variation of SAHE with age 

The mean SAHE and DSAH-lOs are found to have slope 1.92. Figure 6a 
shows the age-dependent variation of mean SAHE and DSAH-lOs for those 
that are computed using all three stations' data. It is seen that both SAHE and as 
such DSAH-lOs bear a cosine function of age. The period is noted as close to 2 
years. A third-order polynomial (spline-fit) of SAHE and DSAH-lOs with age 
could be written as, 

SAHE = 2A3d^-1.6d^+5md + 03 ( 8 ) 

DSAH -lOs = 0.96d^ -2.95d^ +lMd + 0.2 (9) 

A closer look suggests that the reference temperature at the surface and the 
sonde surface temperature used for calculating the surface saturation vapor 
pressure in (5) gives rise to such a sinusoidal variation of SAHE with age of 


sonde. Evaluating the ratio 


esHsr) 

es(Es) 


in (5) by using the Bolton (1980) method and 


using a common reference surface pressure the ratio of vapor pressure 


esr _ exp{(T„a)/(Er+b)} 

ess exp{(7’„a)/(r„ +t)) 


( 10 ) 


or. 


/ ^sr+T„) 
ssHss) 




a(Tsr-Tss) 

(7;, + fc)(r„+fc) 


( 11 ) 


where a = 17.67, and b = 243.5. It is evident from (11) that the ratio of saturation 
vapor pressure in (5) depends upon the differences in measured surface 
temperatures. Figure 7b shows the variation of difference between surface 
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independently measured temperature and the sonde surface temperature with 
sonde age. The nature of variation between the temperature error and SAHE 
variation with age are just opposite (c.f. Figures 7a, and 7b). Hence a decrease in 
the temperature error will increase the SAHE error in RH and vice-versa as in (5). 
If we represent the surface reference temperature and the surface sonde 


temperature as function of age and with a frequency w 

T„=B,Co5(O)rf + 0,) (12) 

T,s=B2Cos((od + (l>2) (13) 

and take the difference in temperature as a cosine function of age 

Trs ACos{(od + 0) (14) 

then, 

Tfs = B2Cos(cod + 02 ) + ACosiokl + 0) (15) 


re-arranging and using the identity Cos[2(.cod + 0^ )] = [icos^ (<urf + ) - 1] and simplifying 
this further for condition 02 = we could write the denominator of (11) as 

(7^^ + h){T^^ +b) = PCos{(od + 03 ) (16) 

hence the right hand side of (11) could be written in terms of the modified 
amplitude and phase 

^^Tsr-Tss) ^ aACos(cod + <l,) ^ p,Cos(wd^4>4) 

+ b){T,, - b) DCosicod -b 03 ) + /^ 

(17) 

whereK = b^ + 

2 

From (17) it is evident that the ratio behaves as a cosine function with sonde age 
but with a combined phase and amplitude. The mean surface temperature error 
could be approximately represented as a function of sonde age using a spline-fit 
of the data as follows: 
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Trs - Ts, = -0.3d^ + 0.9d'^ - 0.7 d + 0.023 


(18) 


d. Improvement in CAPE 

The convective available potential energy (CAPE) was computed using 
the uncorrected soundings and once again after correcting RH data using W02 
algorithm. CAPE is the measure of maximum work done by buoyancy and is 
obtained by integrating the differences in the virtual temperature of parcel (T^. 
parcel) that of surrouudiug (Tv..suit) over all pressure levels in the region 
between the level of free convection (LFC) - where the parcel become positively 
buoyant and the level of neutral buoyancy (LNB) - the height at which the parcel 
regains stability as mentioned in Bohren and Albretch (1998): 

LNB 

CAPE = -Rj J P (19) 

LFC 

The C APE values show a significant increase after W02 based correction is 
applied hus enhancing the site-mean values. Table I gives the quantitative 
information about the CAPE comparison between the old and new sets. 

5. SUMMARY 

We have successfully applied the recently published W02 system-bias 
correction algorithm to all Vaisala soundings obtained during the TRMM field 
campaigns. The day-time and night-time differences in the error (uncorrected 
minus corrected) RH mean profiles for Rebio-Jaru and Rolim de Moura stations 
(LBA sites) are found to be small when compared with similar profiles obtained 
from an oceanic KWAJEX sounding site (R/V RHB). A sonde age (0.06 to 2.06 
years) versus height (up to 650 mb level) plot of corrected minus uncorrected 
specific humidity is shown in Figure 8. It is apparent from Figure 8 that more 
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moisture loading and drying is carried out by W02 algorithm at lower levels 
than at mid-levels and the magnitude of correction increase with age of the 
sonde. However, with sondes aged close to 2 years moist bias correction is 
imparted. Exact reason for this reversal needs to be known. The timely 
publication of W02 algorithm and its simplicity helped enabling us to perform 
dry-bias correction on all the RS80-H data The enhanced CAPE values suggest 
that the soundings represent the LBA and oceanic KWAJ regions in a more 
realistic manner. Humidity corrections might lead to more efficient mesoscale 
analysis in line with the large-scale improvements already observed in the TOGA 
CO ARE data (Ciesielski et al., 2002). However, a surface temperature correction 
is advisable in order to solve for the problems that is found with the 
quantitative SAHE estimation. To check for accuracy of the upper-air data it is 
advisable that the upper-air corrected profiles be compared with some other 
independent measurements taken during the period of TRMM field campaign. 
The version 'e' files with about 367, 462 corrected RH observations as per W02 
algorithm have been posted to Goddard Distributed Active Archive Center (G- 
DAAC). 
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7. FIGURE CAPTIONS 

Figure 1 . The Vaisala sites are marked with an 'X'.a) Locations of Rebio-Jaru 
and Rolim de Moura Vaisala sounding sites shown with respect to 
other LBA VIZ sounding sites in Brazil during January-February, 
1999. (b) Location of R/V RHB Vaisala sounding site during 
KWAJEX at approximately 8.36°N, 67.73®E held during July- 
September, 1999. 

Figure 2 . A comparison plot for the mean profiles of Vaisala measured RH and 
Temperature values and VIZ measured RH and temperature values. 
This plot was obtained after averaging the profile for the entire LBA 
period. Note the reversal of RH bias in the mid-level. 

Fi gure 3 . The system-bias correction algorithm flow-chart based on W02. 

Figure 4 . (a) Profiles of mean daytime and night-time differences between the 

uncorrected and corrected RH data for Rebio Jam (b) Rolim de-Moura 
and (c) RV/RHB site. 

Figure 5 (a) Bar-charts showing observation counts at various sonde age for 
Rebio Jam, site and, (b) for Rolim de Moura site. 
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Figure 6 A Plot showing the variation of GC, CC, SAE, SAKE, and the 

uncorrected minus corrected RH data with sonde age. A sixth order 
polynomial fit is shown for each of the correction parameter involved 
in W02. 

Figure 7 (a) Variation of mean SAHE (and DSAFIlO-s) with sonde age, (b) 

Variation of mean of the difference of surface reference temperature 
and sonde surface temperature with age of sonde. 

Figure 8 Shows the sonde age-height plot of specific humidity corrections 

(g/Kg) imparted to the RS80-H data using W02 algorithm. Note that 
the sign of correction delta changes from negative to positive for the 
sonde age exceeding 1 year. 
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